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Abstract 

We discuss a method of solving n*'* order scalar ordinary differential equations by 
extending the ideas based on the Prelle-Singer (PS) procedure for second order ordi- 
nary differential equations. We also introduce a novel way of generating additional 
integrals of motion from a single integral. We illustrate the theory for both second 
and third order equations with suitable examples. Further, we extend the method to 
two coupled second order equations and apply the theory to two-dimensional Kepler 
problem and deduce the constants of motion including Runge-Lenz integral. 

1 Introduction 

During the past two decades or so immense interest has been shown in identifying in- 
tegrable/solvable nonlinear dynamical systems. Different ideas have been developed or 
employed in isolating such nonlinear ordinary /partial differential equations [1]. Focussing 
our attention on ordinary differential equations (ODEs), one finds that using certain novel 
ideas, Calogero and his coworkers have generated a wide class of integr able/solvable non- 
linear systems and explored their underlying features [2-4]. In this paper we introduce 
another straightforward method of identifying such solvable and/or integrable equations. 
To do so we employ the extended Prelle-Singer (PS) procedure [5]. 

Sometime ago Prelle and Singer [5] have proposed a procedure for solving first order 
ODEs that presents the solution in terms of elementary functions if such a solution exists. 
The attractiveness of the PS method is that if the given system of first order ODEs has a 
solution in terms of elementary functions then the method guarantees that this solution 
will be found. Very recently Duarte et al [6] have modified the technique developed by 
Prelle and Singer [5] and applied it to second order ODEs. Their approach was based on 
the conjecture that if an elementary solution exists for the given second order ODE then 
there exists at least one elementary first integral I{t, x, x) whose derivatives are all rational 
functions of t, x and x. For a class of systems these authors have deduced first integrals 
and in some cases for the first time through their procedure [6]. Recently the present 
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authors have generahzed the theory given in [6] and pointed out a procedure to obtain all 
the integrals of motion/general solution and solved a class of nonlinear equations [7-9]. 

The PS procedure has many attractive features. For a large class of integrable systems, 
this procedure provides the integrals of motion/general solution in a straightforward way. 
In fact this is true for any order. The PS method not only gives the first integrals but 
also the underlying integrating factors. Further, like Lie-symmetry analysis and Noether's 
theorem the PS method can also used to solve linear as well as nonlinear ODEs. In 
addition to the above, the PS procedure is applicable to deal with both Hamiltonian and 
non-Hamiltonian systems. 

In this paper we extend the above theory [6-8,10] to n^^ order scalar ODEs and derive a 
relation which connects integrals of motion with the integrating factors. We demonstrate 
the method to second and third order ODEs. We also introduce a novel way of generating 
all the integrals of motion from a single integral, which is applicable to a class of ODEs, 
and demonstrate our ideas by considering the previous examples. Finally, we prolong the 
method to two coupled second order ODEs. 

We note here that several works have been undertaken earlier to study and classify 
systematically the nth order ODEs based on different ideas. To name a few, besides the 
works of Calogero [2,3], we cite. Lie symmetry approach [11-14], the method of Bluman 
and Anco [15,16], Painleve analysis [1,17] and so on. 

The organization of the material is as follows. In the following section we extend the 
theory of modified PS method to n*'* order ODEs. In Sec. 3, we describe a procedure to 
generate second, third and higher integrals of motion from a single integral, provided the 
first integral can be written in certain specific form. In Sees. 4 and 5, respectively, we 
demonstrate the theory to second and third order ODEs with suitable examples. In Sec. 
6, we briefly discuss the application of PS procedure to higher order ODEs. In Sec. 7 
we prolong the theory to two-coupled second order ODEs and illustrate the theory with 
an example, namely the two-dimensional Kepler problem. We give our final remarks in 
Sec. 8. 



2 Prelle-Singer Procedure 

Let us consider a class of n*'^ order ODEs of the form 

x(") = ^, P,QGC[t,x,xW,x(2),...,x("-i)], (2.1) 

where x*^^) = ^,x^'^'> = ^ and x*^"^ = ^ and P and Q are polynomials in t,x,x^^\x^^\ 
x^"""^-* and x^""^) with coefficients in the field of complex numbers. Let us assume that 
the ODE ()2.1|) admits a first integral I (t , x , x^^\ x^'^\ ...,x^"'~^^) = C, with C constant on 
the solutions, so that the total differential gives 

dl = Itdt + I^dx + I^mdx^^^ + ... + 4(„-i)dx("-^) = 0, (2.2) 

where subscript denotes partial differentiation with respect to that variable. Rewriting 
equation (|2.1I) of the form ^dt — dx^"'~^^ = and adding null terms Si{t, x, x^^^, x*^"~^^) 

x(i)dt-5i(t,x,x(i),...,x("-i))dx and 5i(t,x,x(^), ...,x("-i))xWdt-5i(t,x,x(i), ...,x("-i)) 
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(ix*-* ^\ i = 2,3, ...n — 1, to it we obtain that on the solutions the 1-form 

(j^ + Y^ Six^^ dt - Sidx - S^dx^'-'^'^ - = 0. (2.3) 

i=l ) i=2 



Hence, on the solutions, the 1-forms 1)2.21) and (|2.3() must be proportional. Multiplying 
(tOl by the factor R{t 1)) which acts as the integrating factor for ()2.3p . we 

have on the solutions that 



dI = R 



n—l \ n—1 

+ Y Six^'^ j dt - Sidx - Y S^dx^'-^^ - 



1) 

i=l ■' i=2 



0, (2.4) 



where 4> = P/Q. Comparing equations (|2.2() with ()2.4p we have, on the solutions, the 
relations 



n-l 



Ix — —RSi, 

= -RSi+i, i = l,2,...,n-2, 
4(„-i) = -R. (2.5) 

The compatibility conditions, Itx = Ixt, -^^(n-i) = -^^^(n-i)^, I^xi^-'^) = -^a;("-i)x' ^tx(^) = 
= /a.(i)j,(n-i) = i = 2,3, ... ,n — 2, between the equations 

H2.5|) . provide us the conditions 

(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 
(2.11) 
(2.12) 
(2.13) 

where the total differential operator D is defined by 



D[Si] 


= — 0x + Si(/)^(„-i) + SiSn-1, 




D[Si] 


~ + ^i4'x("-'^) + SiSi^i — S 


'i-i, i = 2,3, . . . ,n - 2 


D[Sn-l] 


= — + 5„_i(/)^(„-i) + S'^_i — 


Sn-2, 


D[R] 


= —R{Sn-l + (/'j.(n-l)), 




R^(i)Si 


= —RSi^i + RxSi^i + iZS'j+i^., ? = 


-- l,2,...,n-2. 




= —RSj^i^(i) + R,j.jSi^i + RSiJ^i^j, 


= 1,2, . . . ,n - 2, 


R^(r) 


= R^(n-i)Si+i + i?S'j+i^{„-i) , i = 1 


2,...,n-2, 


Rx 


= R^(n-l)Sl + RSl^{n-l) , 





9t 9x ^ (9x(*-i)' 

z=2 

We note that (|2.6|1 - (|2.13|1 form an overdetermined system for the unknowns. Si, i = 
1, . . . , n — 1, and R. For example, one can check that for a second order ODE, n = 2, one 
gets three equations for two unknowns, say, S and R, whereas for a third order ODE one 
gets six equations for three unknowns, say. Si, S2 and R. Thus in this procedure, for a 
given n*'^ order ODE, one gets '^^"■^^^ number of equations, for n unknowns, out of which 

"^"2 equations are just extra constraints. 

The crux of the problem lies in solving the determining equations and identifying suf- 
ficient number of integrating factors and null forms. But the point is that any particular 
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solution will suffice for the purpose. We solve equations (|2.(ij) - (|2.13|) in the following way. 
Substituting the expression for ^ = ^, obtained from equation 1)2. into H2.6() - H2.8|) we 
get a system of differential equations for the unknowns Si, i = l,...,n — 1. Solving them 
we can obtain expressions for the null forms Si's. Once Si's are known then equation 1)2. 9|1 
becomes the determining equation for the function R. Solving the latter we can get an 



explicit form for R. Now the functions R and Si 



1, 



1, have to satisfy the extra 



constraints (|2.in|) - (|2.13|) . However, once a compatible solution satisfying all the equations 
have been found then the functions R and Si, i = 1, 2, . . . , n — 1, fix the first integral 



I{t, X, X 



(1) 



,(n-l)^ 



by the relation 



n „ 

/(t,x,x«,...,x("-i)) = ^r,- / 



d 



,i=l 



(2.14) 



where 



rj =-I [RS,-i + ^ (Eill rk)] rfxO-i), j = 3, . . . 



, re. 



Equation (|2.14l) can be derived straightforwardly by integrating the equations (|2.5|) . Now 
substituting the expressions for </>, R and Si,i = 1,2, ... ,n, into (|2.14|) and evaluating the 
integrals one can get the associated integrals of motion. However, we have to point out that 
we have not examined the question of existence of consistent solutions to equations (|2.(i)) - 
1)2.13(1 at present. 



3 Method of generating integrals of motion 

In the above, we derived the integrals of motion, /j, i = 1,2, ... ,n, by constructing 
sufficient number of integrating factors. Interestingly, for a class of equations, one can 
also generate the required number of integrals of motion from one of the integrals, say /i, 
if it is known, provided its form can be written in a specific form. For example, for the 
^th order equation (|2.1() . one can generate I2, 13, ■ ■ ■ , In-i and /„ from Ii itself. In the 
following we illustrate this possibility. 

Let us assume that there exists a first integral for the n*^ order equation (|2.H) of the 
form, Ii{t, 

(" ') = C. Now let us split the functional form of the first 
integral Ii into two terms such that one involves all the variables {t J X J X ^ ^ ^ • • • J X while 
the other excludes x^"~^\ that is, 

h = Flit, X, x^^\...,x^^-'^) + F2{t, X, x«, (3.1) 

Of course, such a splitting is not unique, which can in fact be used profitably further to 
identify new integrals. 

Now let us split the function Fi further in terms of two functions such that Fi itself 
is a function of the product of the two functions, say, a perfect differentiable function 
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■^Gi{t,x,x^^\ ...,x^"' and another function G2{t,x,x^^\ ...,x^" ^)), that is, 

+F2fGi(t,x,x(')....,i("-'))). (3.2 



We note that while rewriting equation in the form H3.2() . we demand the function 

F2(i,x,x(-^\ in automatically to be a function of Gi{t,x,x^^\ ...,x^^~'^'^). 
Now identifying the function Gi as the new dependent variable and the integral of G2 over 
time as the new independent variable, that is, 

w; = Gi(t,x,x«,...,x("-2)), z= rG2(t',x,xW,...,x("-i))(it', (3.3) 



one can rewrite equation ()3.2() of the form 

I = Fi(^^^ +F2{w). (3.4) 
In other words 

Fj^)=I-F2iw). (3.5) 



dz J 

Now rewriting equation 1)3. 5() one obtains a separable equation 

? = /M, (3.6) 

dz 

which can be integrated by quadrature, and I2 can be obtained. 

The procedure given above is easy to follow and can be used to solve a class of problems 
straightforwardly. In fact, for the linearizable second order ODEs, we find that the function 
F2 turns out to be zero and as a consequence one gets ^ = Ii (from (|3.6|) ) which can 
be integrated to obtain the second integration constant trivially. On the other hand, if 
F2 is not zero then the second integration constant can be deduced after the integration 
of (|3.6jl . which can be done for a number of examples. As far as third order ODEs are 
concerned the first order equation 1)3. 6() provides us the second integral whereas a different 
choice Gi and G21 where Gi and G2 are different from Gi and G2, directly leads us to 
the third integral (see Sec. 5). The procedure can in principle be extended to higher order 
ODEs. 



4 Second order ODEs 

The theory developed in Sees. 2 and 3, in principle, can be used to solve a class of equations. 
To illustrate the underlying ideas let us first consider a second order ODE. Since we are 
going to discuss in detail only second and third order ODEs, hereafter, we use the notation 
X, X and 'x instead of x^^\ x^'^^ and x^'^^ for ^ and respectively. 
Fixing, n = 2, the determining equations (|2.6() - H2.13() get simplified to 

St + xS, + cl)Si, = -cl), + Scl)i + S^ (4.1) 
Rt + xR^ + ^Ri = -R{S + (f>i), (4.2) 
R^ = RiS + RSi, (4.3) 
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The integral of motion, (|2.14|) . is fixed by the relation 



I{t,x,x) =ri - r2 

where 

ri 



-R + -7T (n - r2) 

ax 



dx, (4.4) 



= J R{(j) + xS)dt and ^"2 = 1 i^^ + 



The method of solving the determining equations ()4.1() - H4.3|) is described in detail in Ref. 
[6,7]. However, in order to be self-contained, we briefly summarize the main ideas in the 
following. 

Let us first solve the equation (|4.H) with the given (p and obtain expression for S. 
Once S is known then equation 1)4. 2|) becomes the determining equation for the function 
R. Solving the latter one can get an explicit form for R. Now the functions R and S 
have to satisfy an extra constraint, that is, equation ()4.3() . We note at this point that 
all solutions which satisfy equations (|4.1|) - (|4.2j) need not satisfy the constraint (|4.3|) since 
equations ()4.1|) - (|4.3|) constitute an overdetermined system for the unknowns R and S. For 
example, let us consider a set {S,R) be a solution of equations H4.1|) - ()4.2() and not of the 
constraint equation (|4.3j) . After examining several examples we observe that one can make 
the set {S, R) compatible by modifying the form of R as 

R = F{t,x,x)R, (4.5) 

where R satisfies equation (|4.2|) . so that we have 

(Ft + xF, + cl)Fi)R + FD[R] = -FR{S + cj)^). (4.6) 

Further, if is a constant of motion (or a function of it), then the first term on the left 
hand side vanishes and one gets the same equation ()4.2() for R provided F is non-zero. In 
other words, whenever F is a non-zero constant or a function of the integral of motion then 
the solution of equation 1)4. 2|) may provide only a factor of the complete solution R without 
the factor F in equations (|4.5|) . This general form of R with S will form a complete solution 
to the equations (|4.H) - (|4.3|) . In a nutshell we describe the procedure as follows. First we 
determine S and R from equations H4.1I) - (|4.2() . If the set {S,R) satisfies equation (|4.3|) 
then we take it as a compatible solution. On the other hand if it does not satisfy (|4.3jl 
then we assume the modified form R = F{I)R, and find the explicit form of F{I) from 
equation (|4.3|) . which in turn fixes the compatible solution (5, R). Once Ii is derived then 
the second integration constant can be deduced either utilizing our procedure described in 
Sec. 3, or finding another set of solutions {S2,R2) for equations (|4.1j) - (|4.3|) . The method 
has been applied to several interesting nonlinear dynamical systems and interesting results 
have been obtained [7]. In the following example we illustrate both the above ideas. 

4.1 Example 

Let us consider an equation of the following form for illustrative purpose: 



(4.7) 
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so that the equations (|4.H) - (|4.:-{|) become 



(2x-l)±\ _^/x(l+x)-l\.., , 2(2i-l)i„ , „j 



(2x-l)i;2 / 2(2x-iyx\ , , 

R:^ - SRi, - RS± = 0. (4.10) 

As mentioned in Sec. 2, any particular solution satisfying equations H4.8() - H4.1Up is sufficient 
to derive an integral of motion. We solve equations (|4.8j) - (|4.1flj) in the following way. 
Equation 1)4. 8|) is a first order partial differential equation in S with variable coefficients. 
To seek a particular solution for S we consider a simple ansatz for S of the form 

S = a{t,x) +b(t,x)x, (4.11) 

where a and b are arbitrary functions of t and x (in other examples one may need to 
take rational forms in x). Substituting (|4.11j) into (|4.8jl and equating the coefficients of 
different powers of x to zero we get a set of partial differential equations for the variables 
a and b. Solving them we arrive at 



(2x-l) _ 1 , (l-2x ) 

[1 + x^)' ' t^{l + x^ 



S, = )^t—±l±, S2 = -^+'^^^x. (4.12) 



Substituting the forms of Si and 52 separately into ()4.9() and solving the resultant equa- 
tions one can obtain the corresponding forms of R. Let us first consider 5i. Substituting 
the latter into ()4.9|1 we obtain an equation for R of the form 

Rt + iR.+ ^^^^R. =-[^)r. (4.13) 

One can immediately identify a particular solution 
(1 + x^j 

to this equation with a polynamial ansatz in x. One can easily check that Si and Ri 
satisfy equation (|4.1flj) also. As a consequence one can deduce the first integral, using the 
relation 1)4. 4() . of the form 

■ tan~^x 

(1 + x^j 

Now substituting the expression 52 into 1)4. 9|) we obtain an equation for R. As in the 
previous case, one can easily fix a particular solution of the form 

i?2 = --■ (4.16) 

X 

However, this set (52,i?2) does not satisfy the extra constraint (|4.1U|) and so to deduce 
the correct form of R2 we assume that 

R2 = F{h)R2 = -F{h)-, (4.17) 

X 
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where F is an arbitrary function of the first integral Ii. Substituting (|4.17|) into equa- 
tion (|4.1U|1 we obtain F = j^, which fixes the form of R as 

^2 = -7^— (4-18) 

Now one can easily check that this set {S2, R2) is a compatible solution for the set (|4.8jl - 
(|4.10() which in turn provides I2 through the relation H4.4() in the form 

Using the explicit forms of the first integrals Ii and I2, the solution of Eq. (|4.1() can be 
deduced directly as 



x{t) = tan 



log {ht + h 



(4.20) 



However, as shown in the Sec. 3, one can also deduce the second integral from the first 
integral itself. By using the procedure indicated in Sec. 3, by demanding that Ii be put 
in the form 

so that in the transformed variables become w = Gi{t, x) and z = G2{t', x, x)dt' and a 
first order ODE results which can be solved by quadrature. For example, in the present 
case, it is easy to rewrite the first integral ()4.15|1 in the form 1)4. 21() . by inspection, namely, 



d 



h = f^ [^e'^^-'-j , (4.22) 
and identifying with we get 

Gi = e*»""'^, G2 = 1, F2 = 0. (4.23) 
With the above choices, equation 1)3. 3() furnishes the transformation variables, 

w = e*'^'^"'^, z = t. (4.24) 
Substituting into we get 

which in turn gives the free particle equation by differentiation or leads to the solution 
H4.2U() by an integration. On the other hand vanishing of the function F2 in this analysis 
is precisely the condition for the system to be transformed into the free particle equation. 
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5 Third order ODEs 

Now we focus our attention on third order ODEs, ?i = 3 in equation ()2.1() . Fixing n = 3 
in the determining equations (|2.6j) - (|2.13|) . we get 

D[Si] = -^x + SkPx + S1S2, (5.1) 

D[S2] = -(Px + S2(Px + Si - Si, (5.2) 

D[R] =-R{S2 + (l)x), (5.3) 

R'xSi = —RSix + RxS2 + RS2X-, (5.4) 

Rx = RxS2 + RS2X-, (5-5) 

Rx =RxSi + RSix, (5.6) 
where the total differential operator D is defined by 
^ d d d d 
at ox ox ox 
The associated integral of motion is fixed by the relation 

/ = n - r2 - rs - y (^-^ + ^ (^1 ~ - rs)^ dx, (5.7) 

where 

/d 
{RSi + ^^i)*^^! 

rs = JiRS2 + ^{n-r2))dx. 

As mentioned earlier, ()5.1|) - ()5.6() form an overdetermined system for the unknowns. Si, 
S2 and R. We solve the (|5.1|) - (|5.6|) in the following way. Substituting the expression for 
(j) into (|5.H) - (|5.2j) we get a system of differential equations for the unknowns and 5*2. 
Solving them we can obtain expressions for the null forms (5*1, 5*2 ). Once 5*2 is known then 
equation (|5.3)) becomes the determining equation for the function R. Solving the latter we 
can get an explicit form for R. Now the functions R, Si and S2 have to satisfy the extra 
constraints (|5.4() - H5.6|) . 

5.1 Example: 1 

Let us consider an equation discussed by Bluman and Anco from the symmentres point 
of view [16], 

6tx^ 6x2 

x=^^ + —. (5.8) 



i;2 X 



Substituting cj) = ^ + ^ into we get 

Si, + xSix + xSix + + ^) Si, =Si(^^ + m + s2y. 



S2t + xS2x + xS2x +{^ + ^)s2, ='-W + Tr+S2(^+'-f]+S2'-Si, 



X^ X 

6tx^ 6x2 



Rt + xRx + xRx+ {^ + ^]Rx =-R[S2 + ^ + ^ 



(5.9) 
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One can easily verify that equation admits the following solutions, 

S, = 0, 52 = -^^^-t^, R=^, (5.10) 
6tx^ + 4xx 

5i = 0, 62 = , R=^, (5-11) 

x^ 

^ 2x2 etx^ + 2xx ^ ±2 

x^ x^ x^ 

Having determined the functions S'j's and R, i = 1,2, one can proceed to determine 
the associated integrals of motion. Substituting the expressions into 1)5 .Tf) separately and 
evaluating the integrals, one obtains 

h =3t±^ + €, (5.13) 

h =2t±^ + €, (5.14) 

I3 =2x-6t±-€, (5.15) 

respectively. One can easily check that I'^Sji = 1,2,3, are constants on the solutions, that 
is, ^ = 0, i = 1, 2, 3. From the integrals, h, I2 and Is, we can deduce the general solution 
for the equation (|5.8|) straightforwardly. 

In the following we generate second and third integrals, say, H5.14() and H5.15() . from the 
first integral, that is, (|5.13j) . using our procedure described in Sec. 3. 

Rewriting ()5.13|) in the form ()3.2|) we get 

Id, dt dw dw 

xdV ' dz dt dz' ^ ' 

so that 

w = —tx^, z = —X. (5-17) 

Integrating H5.16() and rewriting the latter in terms of the old variables we get 
■ 4 

h = 2t±^ + ^, (5.18) 

X 

which is exactly the same as the one we derived (vide equation (|5.14l) ) earlier through the 
PS procedure. 

To generate I3 from Ii we rewrite the latter in the form (|3.2j) but with different functions 
w and z, namely, 

x^ d , . dt dw dw 

h = —r-nC^x - 3tx) = -— = — , (5.19) 

X dt dz dt dz 

so that 

u; = 2x - 3tx, z = x. (5.20) 
Integrating H5.19() we get 

w = Iiz + h^ h = w - hz. (5.21) 
Substituting (|^:^ and (f^l^ into ^^^T!^ we get 

Is = 2x - 6tx , (5.22) 

X 

which exactly coincides with ()5.15|) . 

In a similar way one can derive Ii and I2 from and Ii and Is from 12- 



V K Chandrasekar, M Senthilvelan and M Lakshmanan 



11 



5.2 Example: 2 

Let us consider another nontrivial example which was discussed by Steeb [18] in the context 
of invertible point transformations, namely, 

X + ^-?,x- — + 2x = {). (5.23) 

X X 

Substituting the form cj) = —^{xx — xx — iP' + into ()5.1|) - ()5.3|) and solving them, we 
find the following three particular solutions which satisfy the determining equations 

=(^,2ii^,xe-2*), (5.24) 
(52,F2,i?2) =(^,^,xe-*), (5.25) 
(53,F3,i?3) =(^^^,^%^,x). (5.26) 

Substituting the expressions (|5.24|) - (|5.26j) separately into the expression (|5.7|1 and evalu- 
ating the integrals we obtain 

h = {x^ + XX - xi;)e~^*, (5.27) 
/2 = {±'^ + XX - 2xx)e-\ (5.28) 
Jg = (^2 + - 3xx + x^), (5.29) 

from which the general solution can be written of the form 

x{t) = (^-^e^' + he' + I^\ (5.30) 
From the integrals /i,/2 and /3, one can also obtain the linearizing transformation. 



6 Extension to Higher order ODEs 

In the previous two sections we discussed the PS procedure applicable for second and third 
order ODEs. Following the same steps one can derive the determining equations for the 
fourth order ODE also from the equations (|2.6() - (|2.13|) . This can be done by restricting 
to n = 4 in these equations. For example, for the present case one obtains ten equations 
for four unknowns, namely, ^i, 5*2, 5*3, and R. Solving them consistently one can obtain 
explicit expressions for the null forms, Sj's, i = 1,2,3, and the integrating factor R. Once 
the SiS and R are known then the associated integral of motion can be constructed using 
the relation 1)2. 14|) with n = 4. Here also, if one is able to find four sets of particular 
solutions, say, {Sij,Ri), i = 1,2,3,4 and j = 1,2,3, one can straightforwardly construct 
four integrals of motion. On the other hand if one has less number of integrals of motion 
then using the procedure described in Sec. 3 one can generate the remaining integrals of 
motion from the known ones and establish the integrability. The extension to higher order 
ODEs follows along similar lines. 
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7 Extension to coupled second order ODEs 

So far we discussed the applicabihty of the PS procedure to scalar differential equations. 
Interstingly, the procedure can also be extended to coupled ODEs. In the following we 
describe the procedure to two coupled second order ODEs and the application to higher 
dimensional equations will be discussed elsewhere. 

Let us consider a class of second order ODEs of the form 

x = ^ = ^, 2/ = ^ = ^> P„Q.eC[t,x,y,x,y], . = 1,2. (7.1) 

Let us suppose that the system ()7.1() admits a first integral of the form I{t, x, y, x,y) = C 
with C constant on the solution, so that the total differential gives 

dl = Itdt + Ixdx + lydy + I±d± + lydij = 0, (7.2) 

where subscript denotes partial differentiation with respect to that variable. Rewriting 
1)7. 1|) in the form 

P P 

-^dt -dx = 0, -^dt -dy = (7.3) 
Qi Q2 

and adding null terms Si{t,x,y,x,y)xdt — Si{t,x,y,x,y)dx and S2{t,x,y,x,y)ydt 
—S2{t,x,y,x,y)dy suitably, we obtain that, on the solutions, the 1-forms 

(— — h Six)dt — Sidx — dx = 0, (7.4a) 
Qi 

+ S2y)dt - S2dy -dy = 0. (7.4b) 

V2 



Hence, on the solutions, the 1-forms 1)7. 2() and ()7.4() must be proportional. Multiplying 
(|7.4a|) by the factor Ri{t,x,y,x,y) and ()7.4b|) by the factor R2{t,x,y,x,y), which act as 
the integrating factors for 1)7. 4ap and H7.4b|l . respectively, we have on the solutions that 

dl = Ri{(pi + Six)dt + i?2(02 + S2y)dt - RiSidx - R2S2dy - Ridx - R2dy = 0, (7.5) 

where 4)i = Pi/Qi, i = 1,2. Comparing equations (|7.5j) and (|7.2|) we have, on the solutions, 
the relations 



It = Ri{^i + Si±)+R2{h + S2y), (7.6a) 

h = -RiSi, (7.6b) 

ly = -R2S2, (7.6c) 

h = -Ri, (7.6d) 

ly = -R2- (7.6e) 
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The compatibility conditions between the equations (|7.(i)) provide us the conditions, 

R2 



D[Si] 

D[S2] 

D[Ri] 
D[R2] 

SlRly 

Rix 

R2y 

Rly 

R2x 

Rly 



2x + -f^Si4'2x + Si(pi± + , 
Ki 



R2 

Ri 

■ -(t>2y - + -^^^'iVly 

-{Ricjjli + R2(t)2x + RlSl), 
: - {R24>2y + Rl4>ly + R2S2)-, 
— RlSly + S2R2X + R2S2X, 

SiRix + RiSix, 

S2R2y + R2S2y, 
S2R2X + R2S2X, 

SlRly + RlSly, 

■ R2x, 



+ s. 



292y 



2> 



dt + ^ix + yiy + 



+ 



where the total differential operator is now defined by 
D 

Integrating equations (|7.6|1 . we obtain the integral of motion, 

d 



I = n + r2 + ra + r4 



R2 + 



dy 



ri + r2 + rs + 



dy, 



(7.7) 

(7.8) 

(7.9) 
(7.10) 
(7.11) 
(7.12) 
(7.13) 
(7.14) 
(7.15) 
(7.16) 



(7.17) 



(7.18) 



where 



ri = J (^Rii(l)i + Si±)+R2{^2 + S2y)Jdt, r2 = -j(^RiSi + ^iri)^dx, 
rg =-j(R2S2 + -liri + r2)]dy, n = - J 



^1 + S ( ^1 + ^2 + r-s 



dx. 



As we did earlier, solving the determining equations (|7.7|) - (|7.1()j) consistently we can 
obtain expressions for the function S'j's and -Rj's, i = 1,2. Substituting them into ()7.18|) 
and evaluating the integrals we can deduce the associated integrals of motion. However, 
unlike the scalar case, the determining equations in the present case are highly coupled 
and pose difficulties to approach them directly. To overcome this problem we adopt the 
following technique. We rewrite equations ()7.7|) - ()7.16() for two variables, namely, Ri and 
i?2) by eliminating Si and 5*2, and solve the resultant equations and obtain expressions for 
Ri and R2 . From the latter we deduce the forms of 5*1 and 5*2 by using the relations 1)7. 9|1 
and (|7.10|) . To implement this algorithm we use all the equations (|7.7j) - (|7.16|) effectively 
such that the functions 5j's i = 1,2, can be eliminated and the functions Ri and R2 can 
be deduced from an optimal set of equations. 

To begin with we deduce the following two identities, 



D[RiSi] =-{Rl(plx+R2(l>2x), 
D[R2S2] =-{Rl(t>ly+R2^2y), 



(7.19) 
(7.20) 
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which can be obtained by combining (|7.7)) - (|7.1()|) . Now we have exphcit forms for the total 
derivatives and i?2'S'2- Let us now take total derivative of equations (|7.9() and (|7.1U|) 

and substitute (|7.19j) and (|7.2()|) in the resultant equations to obtain 

Rm + 2xRux + '^yRlty + '^(plRux + 2(^2^1ty + X^Rlxx + '2'XyRixy + y^Rlyy + (t>ltRlx 
+(l>2tRly + X(j)lxRlx + y(l>lyRlx + '^X(j)iRixx + '^{/(f'lRlyx + 2x(j)2Rlxy + '2y4'2Rlyy 
+x4>2xRly + y4>2yRly + <t>lRlx + 4>2Rly + <t>l4>lxRlx + (plRlxx + 4>2<t>lyRlx 
+<Pl4>2xRly + 4>2<t>2yRly + <t>lx{Rlt + xRix + tjRly + (^iRlx + 4>2Rly) 
+4>lRlyy + 2(t)i(l)2Rlxy + <t>2x{R2t + xR2x + yR2y + <t>lR2x + 4>2R2y) 
-Rl4>lx - R2(t>2x + Rl{4>ltx + X(l)ixx + y(t>lyx + Mlxx + Mlxy) 
+R2{(f>2tx + X4>2xx + y4>2yx + (f>l(t>2xx + 4>2<p2xy) = 0, 

(7.21) 

i?2tt + 2xR2tx + 2yR2ty + 2(^li?2ti + 2(^2-^24?/ + X^ R2xx + 2xyi?2xy + R2yy + (t'ltR2x 
+4>2tR2y + X(t)lxR2x + y01j/-R2i + 2x0li?2ai + 2y(j)iR2yx + 2x4>2R2xy + 2y(f>2R2yy 
+X(t)2xR2y + y4>2yR2y + (plR2x + 4>2R2y + (t>l4>lxR2x + <t>lR2xx + (f>2(plyR2x 
+(t)l4>2xR2y + 4>2(t>2yR2y + (plyiRlt + xRix + yRly + (/'l-Rli + (t>2Rly) 
+(plR2yy + 2(f)i(f)2R2xy + (t>2y{R2t + xR2x + yR2y + (t>lR2x + 4>2R2y) 
-Rl(t>ly - R24>2y + Rl{(t>lty + XCpi^y + {/(plyy + 4'l(plxy + Mlyy) 
+R2{<t>2ty + X(t)2xy + y(t^2yy + (f>l4>2xy + (t>2(t>2yy) = 0. 

(7.22) 

equations (f7T2|l - ([7Tfi|) can also be written of the form 

d d d d 

Rlx = ^(^1*^1)' = -qt{R2S2), Rly = -qt{R2S2), R2x = 'q^^^^^^^' (7-23) 

These identities help us to obtain some additional equations for the variables Ri and i?2- 
For example differentiating (|7.9|) with respect to x and (|7.1()j) with respect to ij and using 
the identities 1)7. 23() in the resulting equations, we get 

Rltx + xRixx + yRlyx + i'lR-lxx + 4'2R2xx + 2Rlx + 2(j)2xR2x + 2(j)i±Ri± 

+R2(t>2xx + Rlhxx = 0, (7.24) 
R2ty + xR2xy + yR2yy + <t^lRlyy + 4'2R2yy + 2i?2{/ + 2(t)2yR2y + 2(piyRiy 

+R24>2yy + Rl<t>lyy = 0. (7.25) 

On the other hand, differentiating 1)7. 9|) with respect to y and using 1)7. 23(1 we get 

Rlty + xRixy + yRlyy + (plR-lxy + 4>2Rlyy + Rly + R2x + 4>2xR2y 

+<t>lxRly + R2(p2xy + Rl4>lxy + (t>lyRlx + 4>2yRly = 0. (7.26) 



The remaining possibility, that is, differentiation of 1)7. 10() with respect to x, leads us to 
the same equation (|7.26|) and so we can discard it. 
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We can further simplify equations (|7.21j) and (|7.22|) by utilizing the equations (|7.24j) - 
fiTI^i . The final form of the equations (f7^ and ((7?^ reads 

Rut + 2xRux + 2yRity + (piRux + hRity + x'^Rixx + 2xyRixy - Ri(pix 



As a result now we have a system of five equations (of course in second order) for the 
unknowns i?i and R2, namely equations (|7.24|) - (|7.28j) . Substituting the expressions for 
<j)i and 02 into (|7.24|) - ()7.28|) and solving them one gets the integrating factors i?iand i?2- 
Once -Rj's, i = 1,2, are known the null forms S'j's, i = 1,2, can be fixed through the 
relation Tn^- TTUi^ . 

7.1 Example: Two-dimensional Kepler problem 

Here we consider a simple, but physically important example, namely, two dimensional 
Kepler problem and illustrate the method developed in the previous section. 
Let us consider the Kepler problem in the x — y plane, that is. 



R2<t>2x + y^Rlyy + (t>ltRlx + <t>2tRly + + V^'lyRlx + XcjyiRixx 

+y(j)lRlyx + Xct)2Rlxy + y4>2Rlyy + X(t)2xRly + y(t>2yRly - (plRlx 
-(t>2R2x + 4>lx{Rlt + xRix + yRly) + Rl{(t>ltx + Xcjyixx + y(plyx) 
+4>2x{R2t + xR2x + yR2y) + R2{4>2tx + Xct)2xx + y4>2yx) = 0, 



(7.27) 



R2tt + 2i:i?2te + '2yR2ty + (plR2tx + 02-^2*1/ + X^ R2xx + 'ixyR2xy - Rl(t>ly 
-R2(p2y + (t>uR2x + 4>2tR2y + ij^ R2yy + X(t)lxR2x + y<PlyR2x + X(j)iR2xx 
+y(t>lR2yx + x4>2R2xy + y4>2R2yy + X(t>2xR2y + y4>2yR2y " 4>lRly 
-<i>2R2y + (f>ly{Rlt + xRlx + ijRly) + Rl{4>lty + X<plxy + y<Plyy) 
+(f>2y{R2t + xR2x + yR2y) + R2{4>2ty + X<p2xy + y(f>2yy) = 0. 



(7.28) 




(7.29) 



where r 



xi + yj and r 



r\. The respective equations of motions are 



X 



Mx,y), 



(7.30) 



X 



(a;2 + ?/2) 2 

y 



4>2{x,y)- 



(7.31) 



y 



{x'^ + y^) 2 



In this case the determining equations (|7.24|) - (|7.28|) simplify to 



Rltx + xRixx + yRlyx + iplRlxx + 4>2R2xx + = 0, 

R2ty + xR2xy + yR'2yy + ^l-Rlyy + 4'2R2yy + 2i?2j/ = 0, 
Rlty + xRixy + yRlyvy + 4'lR'lxy + <p2Rlyy + Rly + R2x = 0, 



(7.32) 

(7.33) 
(7.34) 
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Rltt + 2xRux + '^yRlty + X^Rlxx + '^XyRixy + y^Rlyy + i<Ala:-Rli + {/(plyRlx 
+X(t)lRlxx + y(t>lRlyx + X(t)2Rlxy + y(t>2Rlyy + X(l)2xRly + y(t>2yRly 

+4'lRltx + 4'2Rlty - Rl(t>lx - R24>2x - 4>lRlx - 4>2R2x = 0, (7.35) 
R2tt + 2xR2tx + '2yR2ty + X^R2xx + '2xyR2xy + y^R2yy + X(pixR2x + y(plyR2x 
+X(t)iR2xx + y(t>lR2yx + X<p2R2xy + ^02-^2^ + X(l)2xR2y + y(l>2yR2y 

+4>lR2tx + (t>2R2ty - Rl4>ly - R24>2y - 4>lRly - 4>2R2y = 0. (7.36) 

To solve equations 1)7. 32() - (|7. 36(1 we seek an ansatz 

Ri = ai{x,y) + a2{x,y)x + a3{x,y)y, (7.37) 
R2 = bi{x, y) + 62(2;, y)x + 63(x, y)y, (7.38) 

where Oj's and 6j's, i = 1,2, 3, are arbitrary functions of x and y. Substituting ()7.37|) and 
1)7. 38(1 into H7.32() - H7.36() and solving the resultant equations we can obtain at least the 
following three solutions. 

(i) Ri=x, R2 = y, (7.39) 

(ii) Ri = y, R2 = x, (7.40) 

(Hi) Ri = 2yx — xy, R2 = —xx. (7.41) 

We are also now searching for other possible forms with a modified ansatz. Substituting 
TTMi into (IZini) and (tTTTHl we get 

(z) Si = ^, 82 = ^ ^. (7.42) 

i;(x2 + y2)2 y(x2 + 2/^)2 

In similar way equations (|7.4Uj) and (|7.41l) provide us 

{ii) Si = -y-, 52 = --, (7.43) 

y X 

(xy \ /' ' 2 \ 

-xy+-pj) ~ (X - pr) /^7/l/l^ 

[in) ^1 = — 52 = r^. (7.44) 

2yx — xy xx 

Substituting the expressions i?i's and 5j's, i = 1,2, into (|7.18j) and evaluating the integrals 
we get, 

I, = h±^ + f)- (7.45) 

which is of course the Hamiltonian of the system. The forms -Rj's and 5i's, i = 1,2, provide 
us 

h = yx- xy, (7.46) 

the second integral, namely, the angular momentum. The integrating factors -Rj's and null 
form 5i's, i = 1,2, lead us to 

h = x{yx-xy) , (7.47) 

namely, the Runge-Lenz constant. 
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8 Final Remarks 

In this paper, we have discussed the method of solving a class of ODEs through the mod- 
ified PS method. The method is applicable to both scalar and multicomponent equations 
of any order. We also demonstrated the theory with examples. Apart from the above, in 
the scalar case, we introduced a novel way of generating integral of motion from a single 
integral and illustrated our ideas with the same example considered previously. The appli- 
cation of this method to multicomponent systems and their integrability and linearization 
properties will be published elsewhere. 
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